function [c,r] = findbrightest(I)
    close all;
    Iorg = I;
    [sx,sy,~] = size(I);
    h1 = fspecial('gaussian',[sx sy],(sx+sy)/10);
    
    hsv = rgb2hsv(I);
    v = hsv(:,:,3);
    
    h = fspecial('gaussian',10,5);
    v = imfilter(v,h);
    
    g = gradient(double(v));
    
    g = g .* h1;
    
    figure; imshow(g)
    max_g = max(g(:));
    
    [r, c] = find(g == max_g);
    
    imshow(Iorg);
    figure; imshow(Iorg); hold on;
    plot(c,r,'ro','MarkerSize',10);
end